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Abstract 

We take account of the lowest order back reaction on the fluid and 
modify the Kulsrud and Anderson equation dtS-u = ^l^M obtained 
in the kinetic dynamo theory, where Em is the energy density of the 
magnetic field. Furthermore, we apply our present result to some as- 
trophysical stages where the magnetic field is expected to be amplified 
by the dynamo mechanism. 
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1 Introduction 



The magnetic fields have been observed in various astrophysical scales [|l|]. 
The origin is one of the important problems in cosmology ||. Although 
an attractive mechanism in protogalaxy was proposed by Kulsrud et al@ 
for the galactic magnetic field, it cannot explain how the magnetic field is 
made in intergalactic and intercluster regions^. Thus, it is worth inves- 
tigating the generation and the evolution of the primordial magnetic field. 
Since the strength of these fields strength is too small we expect that the 
amplification of these fields occurs due to the dynamo mechanism. It is well 
known that the mean magnetic field can be amplified enough to explain the 
present observation in the kinetic dynamo theory Q. However, as Kulsrud 
and Anderson showed H, the growth rate of the fluctuation around the mean 
magnetic fields is much larger than that of the mean field in interstellar medi- 
ums. This means that the kinetic dynamo theory breaks down. Hence, one 
must investigate the effect of the back-reaction to the kinetic theory. So far 
the back-reaction on the mean field has been considered §. Setting apart the 
problem of the kinetic dynamo theory in interstellar mediums, it is obvious 
that the kinetic theory cannot hold near the equipartition state in general. 

In this paper, we consider the back-reaction on the fluctuation and derive 
the evolutional equation of energy of the magnetic field, that is, modified 
Kulsrud and Anderson equation. Then we apply it to some examples. 

The rest of the paper is organized as follows. In Sec. II, we derive 
the modified Kulsrud and Anderson equation with the lowest order back- 
reaction under a phenomenological assumption. In Sec. Ill and IV, we give 
applications and remarks, respectively. 

2 Modified Kulsrud and Anderson Equation 

In the Fourier space the basic equations of the incompressible MHD are 

— — |^(k- q,t)v*(q,t) -6j-(k-q,t)6fc(q,t)J (1) 



and 



dA(k, t) = ikj J [«i(k - q, t)bj(q, t) - Vj(k - q, t)&i(q, £)] , (2) 
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where bi(k, t) := Bi(k,t)/ y/4np, p is the energy density of the fluid and 
Pjjfc(k) = kjPik(k) = kj(5ik — /cj/cfc/|k 2 |). For simplicity, we neglected the 
diffusion terms in the above equations. This simplification can be justified 
by the fact that almost of astrophysical systems have high magnetic Reynolds 
number. 

Following Kulsrud and Anderson |J and considering a small time step as 
the parameter of expansion, we evaluate the time evolution of the magnetic 
field by iterations; 

k(k,t) = b t (k,0) + bf\k,t) + b ( ?\k,t) + --- 

= fof ) (k) + 6? ) (k,t)+6f ) (k, *) + •--, (3) 

where b\ ^ (k) is the initial field. For the fluid velocity, we take account of the 
back-reaction from the magnetic field (Lorentz force) as follows; 

v l (k,t) = vl 1 \k,t)+6v t (k,t), (4) 

where i>| (k, t) is statistically homogeneous and isotropic component and 
satisfy 

(^teOuf^q,*)) = (2nf[j 1 (k)P tJ (k)+ i J 2 (k)e lkj k k ]6 3 (k-q)6(t-t') 

= (2n) 3 V lJ (k)5 3 (k-q)5(t-t'). (5) 

This statics holds in the region where is far from the boundary. J\{k) and 
J2(k) denote the velocity dispersion and the mean helicity of the fluid, 

(v«(x,t)-v( 1 )(x,t)) = 2| J^Lj^kMO) 

and 

(v«(x,f) • V x v( 1 )(x,t)) = -2 J ^Lk 2 J 2 (k)5(0). 

The second term in the right-hand side of the eq.(4), 5vi(k,t), is determined 
by eq. (1) and this corresponds to the back-reaction term from the magnetic 
field. 

In each orders, the MHD equation becomes 

d t b?(k,t)^2ik j J ^ % (k-q,t)^(q) (6) 
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d t bf\k,t) = / ^^(k - q, (7) 



(27T) 

and 

^(M) = 2iP ijfe (k) / ^^(k - qOftg^q,*)- (8) 

The last equation contains an effect of the lowest order back-reaction on the 
fluid and it gives an explicit expression of Svi(k,t) 

5 Vi (k,t)^2iP ijk (k) [dt'J ^ 6 g>( k -q)&g>(q,0. (9) 

From the eqs. (6) ~ (9), the time derivative of the energy becomes 
<9 t (|b(k,t)| 2 > = (6l 1) *(k,t)^ ) (k,t))+6f ) *(6p ) (k,t)) + c.c. 

= 4 [dt'J ^^fc^*(q)(T,3(k-q,O«[i(k-p,0>6S ) (p) 

-4 [dt f ^^fc i ?€6i° ) ^k)(t;f i (q-k,0T;[,l(q-P,< , )>^ ) (p) 
+ c.c, (10) 

where 

«(k,i'H(q,0> ^ (^(k.O^Cq,*)) 

+(t;{ 1) *(k,0^-(q,*)> + <^*(k,t>f( q ,t)> 
= (2vr)%(k) ( 5 3 (k-q)5(t-t / ) 

~ 4 jf (03 P ^^M2(p - k)^](k)6^(q- P) 

-4 /V/ ^^M^p-^^w^^^-p) 

=: (27r) 3 Vy(k)<f 3 (k - q)5(t - f ) + <5«(k, t>,(q, *)>■ (") 

The above eq. (10) with (11) is the formal equation with the effect of the 
back-reaction. 

Let us consider the simple example with the following initial condition 
6l 0) (x)=M« or bf ) \k)=b (2n) 3 5 3 (k)5 lz . (12) 
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This condition holds approximately as long as the spatial scale of the mag- 
netic field is much larger than the typical scale of eddies. In this case, the 
eq. (10) becomes 

d t (\b(k,t)\ 2 ) = 2(27c) 3 5(0)k 2 z V lt (k) + 2 J* dt'k 2 z 5(v*(k,t')v t (k,t))b 2 

= A{27:fklJ l {k)bl5\Q) -6(2n)XbU^) 2 Mk)5%0)(l3) 

where (A£)& is the time scale of the eddy turnover and its expression will 
be given below. We assumed that the time integral should be estimated as 
J Q d1/[- ■ ■] ~ (At)fc[- • •] in the second line of the right-hand side of the eq. 
(13) because the back-reaction works only during the time scale of the eddy 
turnover. Here we assume Kolmogoroff spectrum for the inertial range f\ 
k < k < fc max ~ i? 3//4 /c o 0, k is the wave number of the largest eddy and R 
is the Reynolds number. From the definition of the velocity dispersion 

/fj Jc Z'&max 
— 7^)5(0) =: dkl(k), (14) 

we obtain the relation 

Hk) = ^JiCWAfk 1 * \v 2 ^ (15) 

where vq is the typical velocity(t>o ~ \J (v 2 )) and we used 5(0) ~ (At)^ 1 . The 
expression of (At)k is given by the estimation of the order of the magnitude 
in the eq. (14), that is, (l/k(At) k ) 2 ~ kl(k). 

Integrating the above equation (13) over k, we obtain the modified Kul- 
srud and Anderson equation 

d t pM = 27p M - 2Cp M , (16) 

where 

2 The inertial range is denned by the scale which is smaller than the largest eddy (~ k^ 1 ) 
and larger than a small scale (~ R^ 3 ^ 4 ^ 1 ) under where the viscosity term is dominant. 
In this range, the transfer of the energy works from large eddy to small one without 
the dissipation of the energy. This leads a sort of 'equilibrium state' with Kolmogoroff 
spectrum^] (Kolmogoroff Theory). 
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7 := 2 / (0s*^i(*) (18) 



and 

t-.= zj^k%{k){M)l (19) 

In the above derivation, we used <5 3 (0) ~ V, where V is the typical volume 
of the system. 

Now we evaluate the coefficients 7 and (. Results are given by 
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dkk 2 I(k)(At) k ~ / ^ dk[kl{k)} 1/2 ~ vok^k 2 ^ ~ i? 1/2 v £<) (20) 

feo Jka 



and 

C~ / dkk 4 I(k)(At)l~ / <i/c[/c/(£;)] / ~ ~ -R— , (21) 

respectively. Defining a dimensionless quantity fiM '■= Pm/vq, we can see 
that the eq.(16) becomes 

d t /i M = 2j/i m - 2(74, (22) 

where (' = (vq ~ Rk^v^. The second term in the right-hand side of the eq. 
(22) comes from the effect of the back- reaction effect. One can see easily 
from the above equation that the back-reaction gives an opposite effect to 
the original kinetic term and make the energy of the magnetic field balance 
with the energy of the fluid. 

Although we know from the procedure used here that the eq. (22) holds 
only in a small time step as pu < 1, we try to extrapolate. As a result we 
find the solution 

m = ?i_(i_jL:rW (23) 

One can see easily that the magnetic 'energy' goes toward the terminal value 
fx* M = 7/C ~ i? _1//2 for a time scale ~ 7 -1 . This value corresponds to the 
saturation value, which is estimated naively on the assumption that the drain 
by the magnetic field is comparable to the turbulent power[]3|]||. 
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3 Applications 



In this section we apply the eq. (22) to two examples which the magnetic field 
is amplified by the dynamo mechanism. First, we treat the time evolution 
of the magnetic field during the first order phase transition in the very early 
universe. We also consider briefly the amplification of the magnetic field in 
interstellar mediums. 

3.1 Electroweak Plasma 

There are attractive mechanisms of the generation of the primordial magnetic 
field in the course of cosmological phase transitions ||. In these scenarios the 
strong magnetic field is expected to be amplified by MHD turbulence during 
the first order phase transition. The detail of the amplification has been 
discussed by using Kulsrud and Anderson equation in the ref. |§. 

We reconsider the amplification of the magnetic field during the phase 
transition by using the modified Kulsrud and Anderson equation. The time 
scale for the equipartition is t equ i ~ 7 _1 ~ R~ x ' 2 VQ X k$ 1 . Since the Reynolds 
number is R ~ 10 2 |J, we can see that it is the same order with the time scale 
of the phase transition. Thus, the magnetic field can be amplified enough 
and the final energy is given by 

E* M ~ R~ 1/2 S V ~ 0.1 x S v , (24) 

where £ v is the energy of the plasma fluid. 

3.2 Interstellar Mediums 

As we stated in Introduction, the kinetic dynamo theory breaks down in in- 
terstellar mediums ||. For interstellar mediums, typical values of key quan- 
tities are 2n/k ~ lOOpc, t> ~ 10 6 cm/s and R ~ v^jk^v ~ 10 8 , where v 
denotes the kinetic ion viscosity; v ~ 10 18 cm 2 s _1 @. Then the typical time 
scale is given by tisM ~ 7 _1 ~ 10 2 yr. Since the time scale of the mean field 
is ~ 10 10 yr[]5]], we realize again the mean field theory is meaningless in the 
present perturbative approach. The final energy of the magnetic field is given 
by^~l(T 4 x£ v . 
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4 Concluding Remark 



In this paper, we considered the lowest order back reaction to the kinetic 
dynamo theory and modified the equation for the energy of the magnetic 
field. As a result we obtained the successful time evolution of the energy 
of the magnetic field. That is, the terminal value of the magnetic energy 
obtained from the eq. (23) equals to the previous qualitative estimation of 
the saturation energy 0||. We also presented the expression depending on k 
(eq. (13)), with which we can evaluate the evolution of the magnetic field for 
various scales. Since the present formalism is general, our equation is useful 
for other situations, for example, the fireball model for 7-ray bursts [TO]. 

Finally, we should comment on our assumption for the initial condition 
(eq. (12)) and the extrapolation of the eq. (22). We choose the initial 
condition in order to obtain the simple result like the eq. (22). Although 
this assumption holds approximately in some cases, it may not be correct 
in general cases. We should also note that we considered only the effect 
of lowest order back-reaction. Properly speaking, if one wishes to analyse 
the vicinity of the equipartition, one must take account of effects of higher 
order back-reaction. The study near the equipartition might become clear 
by using something like the renormalization group approach. The study for 
more general initial condition and with higher order back reaction should be 
done in the future. At the same time, the spatial structure as the typical 
coherent length of the magnetic field also should be discussed. 
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